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Abstract 

We consider the correlation functions of the tachyon vertex operator of the super 
Liouville theory coupled to matter fields in the super Coulomb gas formulation, on world 
sheets with spherical topology. After integrating over the zero mode and assuming that the 
s parameter takes an integer value, we subsequently continue it to an arbitrary real number 
and compute the correlators in a closed form. We also included an arbitrary number of 
screening charges and, as a result, after renormalizing them, as well as the external legs 
and the cosmological constant, the form of the final amplitudes do not modify. The result 
is remarkably parallel to the bosonic case. For completeness, we discussed the calculation 
of bosonic correlators including arbitrary screening charges. 
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1- Introduction 



Two dimensional gravity is not only a toy model for the theory of gravitation, but also 
describes phenomena such as random surfaces and string theory away from criticality 1 . 
The discretized counterpart, namely matrix models, proved to be an efficient means to 
obtain information, especially about non critical string theory, while computing correlation 
functions 2 . 

In the continuum approach 4 , in the conformal gauge, we have to face Liouville theory 5 . 
However, although several important developments have been achieved 6 , we still lack some 
important points, in spite of much effort which has been spent. 

In particular, it is difficult to calculate correlation functions in a reliable way because 
perturbation theory does not apply. Recently, however, several authors 7-14 succeeded in 
taming the difficulties of the Liouville theory and computed exactly correlation functions 
in the continuum approach to conformal fields coupled to two dimensional gravity. The 
technique is based on the integration over the zero mode of the Liouville field. The resulting 
amplitude is a function of a parameter s which depends on the central charge and on the 
external momenta. The amplitudes can be computed when the above parameter is a non- 
negative integer. Later on, one analytically continues that parameter to real (or complex) 
values. The results for the correlation functions of the tachyon operator thus obtained 
agree with the matrix model approach. 2 General correlation functions including arbitrary 
screening charges are also computed, and may be useful for the purpose of studying fusion 
rules of minimal models coupled to two dimensional gravity. 25 

Opposite to the bosonic case super symmetric matrix models are up to now scarcely 
known. There are few papers in the literature concerning this approach 3 ' 25 , and it is desir- 
able to have results from an alternative approach for comparison. For this reason, several 
groups are studying the continuum approach to two dimensional supergravity 17-21 . Our 
aim here is to investigate the supersymmetric Liouville theory. We shall compute supersym- 
metric correlation functions on world sheets with spherical topology in the Neveu-Schwarz 
sector, where the super-Liouville is coupled to superconformal matter with central charge 
c < 1, represented as a super Coulomb gas 23 ' 17 . The results are remarkable, and very paral- 
lel to the bosonic case; after a redefinition of the cosmo logical constant, and of the primary 
superfields, the resulting amplitudes have the same form as those of the bosonic theory 
obtained by Di Francesco and Kutasov 9 . Our present results generalize those presented in 
a recent paper 19 , as well as others recently obtained in the literature 20 ' 21 ' 26 . 

The paper is organized as follows. In section 2 we review some computations of the 
bosonic correlators and generalize them to include an arbitrary number of s.c. . In section 
3 we calculate the iV-point correlator of the Neveu-Schwarz vertex operator in the n = 1 
two dimensional supergravity coupled to N = 1 supermatter including an arbitrary number 
s.c. 

Our results include the limiting cases c = 1 (c = 1) and in those situations some 
physical conclusions can be drawn in the section 4: the amplitude factorizes, and the 
expected intermediate poles have zero residue, due to strong kinematic constraints. All 
results are at least consistent with the matrix model approach, with possible exception 
about the inclusion of screening charges, in the 3-point correlation function. In the latter 
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case, fusion rules deserve careful study. 



2- Bosonic Correlation Functions 



Some iV-point tachyon correlation function in Liouville theory coupled to c < 1 con- 
formal matter were recently calculated by Di Francesco and Kutasov 9 . They worked in 
the DDK's framework 22 where the total action is given by: 



g ab d a( pd b( p - + 2^ + g ab d a Xd b X + ^RX 



(2.1) 



here <p represents the Liouville mode and X is the matter field with the central charge 
given by c = 1 — 12a 2 ). From the literature 17 we know that the constant Q is determined 
by imposing a vanishing total central charge, and is given by 



Q = 2^2 + a%. 



The value of a is determined by requiring e Q< ^ to be a (1,1) conformal operator, yielding 
the equation ^(a + Q) = 1, whose solutions are labeled by a± 

a± = — — ± | cto | i ct+a- = 2 (2.2) 

the semiclassical limit (c — > — oo) fixes a = a+. 

The gravitationally-dressed tachyon amplitudes are the objects we are interested in: 

(T kl ■ --T kN ) = (j[J rf 2^. e ^x( z ,) +/ 3(fe,)^( z ,)^ (2 3) 

where we fix the dressing parameter (3 imposing e lk j x +l3j<l> to be a (1,1) conformal operator 
and supposing the space-time energy to be positive: 

E = p(k) + Q = \k j -a \ . (2.4) 

In the calculation of the amplitudes (T kl ■ • • T kN ) the main ingredient is the integration 
over the matter (Xq) and the Liouville (<po) zero modes. This is the so called zero mode 
technique: one splits 4 ' 5 both, the matter and the Liouville fields as a sum of the zero 
mode (^o), (</>o) plus fluctuations (X), ((f)), where the fluctuations are orthogonal to the 
zero mode. After such splitting we are left with the following integrals: 

VX e lXt) C-i fc *- 2a °) = 2nd ( ^ h - 2a ) , 
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where we have used that on the sphere ^ J d 2 w\fgR = 1 and 




(2.6) 



We thus obtain for the amplitude 

(T kl ■ ■ ■ T kN ) = 2*6 ( k 3 - 2a ) ^ (£) 



N 

X 



(2.7) 



3 = 1 J ^ 



where (• • -)o means that now the correlation functions are calculated as in the free theory 
{ji = 0). The strategy to obtain An is to assume first that s is a non-negative integer and 
to continue the result to any real s at the end. Thus, using free propagators 

(X(w)X(z)) = ((f>{w)(f>(z))o = ln\w- z\~ 2 (2.8) 

and fixing the residual SL(2, <D) invariance of the conformal gauge on the sphere by choosing 
(zi = , Z2 = 1 , zs = oo), the 3-point function is written as: 

Mh,k 2l k 3 ) = ^^^y [ f[<Pw j \w j \ 2a \l-w j \ 2 ef[\w i -w j \ 4p , (2.9) 

a+ 71 j=l i<j 

where we have defined a = —a + (3i , (3 = —a+fo , p = —a+/2. Choosing the kinematics 
fei, ks > cko , k,2 < OiQ < (notice that our notation differs from Ref.[9] by the exchange 
of ki and k^.) we can eliminate (3 using (2.4), (2.6) and the momentum conservation, one 
can write the 3-point amplitude in a rather compact form 

A, = [^{-p)] s \{/\(\{(3 2 -k 2 )) , (2.10) 

j=l V / 

where A(x) = F(x)/T(l — x). After redefinitions of the cosmological constant and of the 
external fields as 

^Ari) ' Tk ^ ^m-k 2 )) ' (2 - n) 

Di Francesco and Kutasov 9 obtained for the three-point function 



A 3 = H' 

5 



(2.12) 



which is also obtained in the matrix model approach. In the following we shall see that a 
similar expression holds for general iV-point tachyon amplitudes with an arbitrary number 
of screening charges. 

The screening charges are introduced in the form of n operators e ld+x and m operators 
e ld - x , with d± solutions of: \d{d — 2a ) = 1 , (d + d- = — a + ct- = —2). Integrating over 
the zero-modes again we get: 



(2.13) 



«vww fen/***"") fen/*'*™ 

3 

= 2n5(J2 h + nd + + md - ~ 2« )^3 m (^i, k 2 , fc 3 ) 
i=i 

where the amplitude A^ m (ki, k 2 ,k 3 ) is given by the expression 

a+ n i=l J i<j 

m „ m 

n J dWn^'ll-nl 2 ?' U ln-r^' (2.14) 



m ^ m 

X 

i=l ° i<j 
n m s 

X 

i=l j=l i=l " i<j 



7T Tin s s 

nn ' ^r 4 n / - n i* - * 



.|4p 



The parameters a, (3 and p are defined as before, and the remaining parameters are 

1 2 

a = d+ki , j3 = d + k 2 , p = -d + 

1 (2 ' 15) 
a' = d-ki , /?' = <i_/c2 , p' = ttgP 

Notice that the gravitational part of the amplitude (integrals over zi) is the same as 
in the case without screening charges. The integrals over U and rj (matter contributions) 
have been calculated by Dotsenko and Fatteev 16 (See their formula (B.10) in the second 
paper of [16]); the result turns out to be 

m n 

A nm = ^ S r(-s)r(s + l)n s+n+m p-* nm [A(l-p)] n [A(l-p'] m JJ A(ip'-ra) Yl A(ip) 

i=l i=l 

m — 1 

x J'[ A(l-n + a' + ip')A(l - n + ft + ip')A(-l + n - a' - ft - (n - 1 + i)p') 

i=0 
n-1 

x Y[ A (l + ol + ip)A(l + (3 + ip)A(-l + 2m-a-j3-(n-l + i)p) 



i=0 

s s — 1 

i=l i=0 



- p)] s n A(ip) Y[ A(l + a + zp)A(l + /3 + ip)A(-l-a-^-(s-l+i)p) 

(2.16) 
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Assuming ct < 0, and the same kinematics (ki , ks > cto , k 2 < cto) we can eliminate 
(3 again, using (2.4), (2.6) and momentum conservation. In this way, we have 

a = a — 2p , a! = —2 + pa 
13 = — 1 — m — (s + n)p , 

(2-17) 

(3 = m — 1 + (s + n)p , (3 = s + n + p (m — 1) 
p=-p , p' = -p~ x . 

Substituting in (2.16) we obtain the three point function with arbitrary screenings 

A nm = ^yr(-s)T(s + l)Tr s + n+m (p)- 4nm [A(l + p' 1 )} ™ [A(l + p)] n 
m n s 

n a^- 1 ~ n ) n a(-<p) [A(i - p)r n am 



m 

X 

i=l i=l i=l 

n— 1 s— 1 

X 

i=0 i=0 
m—1 



| J A(to + (s + n — i)p) Y\ A(— m — (s + n — i)p) 

i=0 i=0 
m — 1 

x ] [ A(i + s + (m - 1 - Op -1 ) 

i=0 

71-1 

J [ A(-l - n + /9 _1 a - z/9 _1 )A(l - s - p _1 a + ip~ x ) 

i=Q 
i-l 

J\ A(l + a - (i + 2)p)A(to - a - (s - 1 - i)p) 

i=0 
s-1 

] J A(l + a + ip)A(m - a + (n + 1 - i)p) . 



i=0 
m — 1 

X 

i=0 
n-1 

X 

i=0 
s-1 

X 
i=0 



(2.18) 



To get a simpler expression for A^ m we look for the term A(p-a) A(p(s^+l)+a-^rn+l) x 
A(-m/9 _1 -(s+n)) which corresponds to Yl i=1 A(|(/??— fcf )). After algebraic manipulations 
we get 

= [^A(-p)] s [-ttA^" 1 )] 7 " [-^A(p)] n f[(-7r)A Q(/3 2 - kf)) . (2.19) 

i=i ^ ' 

This result has been also obtained by Di Francesco and Kutasov 9 ' 20 , as well as Aoki 
and D'Hoker 14 . Note that the factors A(p) and A(p _1 ) can be easily understood; the 
screening operators are renormalized like the tachyon vertex operators with vanishing 
dressing (3(k). 

Remembering the momentum conservation law: 

3 

= (1 - + (1 - m)d_ , (2.20) 



i=l 



we can use (2.19) with n = m = 1 and the formula below to obtain the partition function 
Z: 

°" Z =Al 1 (k i ^0) (2.21) 



d/j 3 
thus, 

We can also get the two-point function by taking, e.g., k 3 — > in (2.19) and (2.21): 

^r(fci,fc 2 ^3^o) = |-^r(fc 1 ,fc 2 ) , (2.23) 

where (n, m) are fixed by (2.20) with k% = 0. Thus we finally arrive at 

p _1 (^Ei=iA + p-iJ i=1 v z / 

(2.24) 

We are able now to calculate ratios of correlation functions to compare with other 
results in the literature. We have (in a generic kinematic region): 

R= A^ m (k u k 2 k 3 )Z Ui=i a +\ k i- a o\ (225) 

A^ikukJAT^ik^k^A^^h) (p-l)(l + p) 

In the case of minimal models the momenta assume the following form: 

k ri< = ^^d+ + ^-^-d- . (2.26) 
Plugging back in (2.25) we find for the ratio R the result 

{T rir 'T r2r 'T r:}r ' 3 ) 2 Z f a +\ 3 nLi \ r i d + + r i d - 



R 



(T riri T rir /)(T r2r /T r2 ^)(T r3 ^T r3 ^) V 2 / (p-l)(p+l) 



(®+Y lU=i \na+ + r l a-\ 
UJ (p-1)(p + 1) {2 - 27) 



That is exactly the same result obtained by Dotsenko 10 (in a given kinematic region) 
which reduces to the result of Goulian and Li for r« = r[; all these results agree with those 
obtained by the matrix model approach. 

At this point the following remark is in order. After renormalizing the screening 
charges: 

p id-X p id-X 
id+X f id^X (2 2%) 

A(p) ' A(p-i) ' (2 ' 28) 

and the tachyon operators as well as the cosmo logical constant \i as before (see (2.11)) 
we get the renormalized amplitude. 

AT = /I s . (2.29) 



Using the above result we would be able to exactly reproduce the ratio in (2.27). 
Thus, the comparison of those ratios with matrix models is not a precise test and only 
measures the scalling of the amplitudes w.r.t. the cosmological constant. All singularities 
contained in the A functions cancel out in such ratios. If on the other hand compare, 
the 3-point function (2.19) directly with the matrix model result we would find a precise 
agreement only for c = 1 («o = 0) where the amplitude can be written as a function of the 
renormalized cosmological constant as: 

3 

A% m ~W]\rO--y/2\ki\) • (2-30) 

i=l 

For c < 1 the amplitude obtained via matrix model are finite and the singularities 
contained in (2.19) are not observed. Neither do the fusion rules for minimal models (see 
discussion in [27,28]) appear in (2.19), although they can be seen via matrix models. 

We now generalize the cases known in the literature for N > 3 including screen- 
ing charges. Repeating the zero-mode technique in the most general case of an iV-point 
function with arbitrary screening charges we have 



AT = (-)"(£) > • 11/ II / ,' 11 



i=l" 3 = 1" k=l 



ml 

(2.31) 

X Yl I d 2 Wi ^iX("i) e id+X(ti) e id- x (r n )^ ^ e 0i*(*i) e «+*(«>i)' 
1=1 J ° 

where s = — -^-(J2iLi A + Q) an d a factor n 3 /a + has been absorbed in the measure of 
the path integral. Fixing the SL{2, (D) symmetry we get: 

rnm 



AT = (-nf {^) S T(-s)I 

N 

J A A ^ A — 



\ 7 . _ 7 .\ A Pij 
\*t *>3\ 



j=4 i<j=4 

/s s s N 

nrfVKi 2a ii-^i 2/3 n \wi-wj\ 4p iii[\wi-z3\ 2p3 

i=l i<j=l i=lj=4 

n n n N (2.32) 

x n d%\uni - uf? n \u ~ t 3 ft n n \ Zj - ^ 

i=l i<j i=l j=4 

m N 

2a' 



/lit iil lit i y 

n^'v 2 "'' n ' v " n n ' 
i=l i<j=l i=l j=4 



3 



n m 



< nn ' 

i=ij=i 
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where 

ctj = kikj — PiPj , o.j = d+kj 

Pj = k 2 kj - fcPj , OLj = d-kj (2.33) 

Pij = ^(hkj - Pi0j) ,Pj = -a+Pj , 4 < j,l < N . 

The integral above has been calculated by Di Francesco and Kutasov 9 , for the case 
n = m = 0. For arbitrary n, m we shall use the same technique. Notice that translation 
invariance 

Wi -> 1 - Wi , Zi -> 1 - Zi , -> 1 - t, , r 2 -> 1 - r, 
implies the symmetry relations 

a <-> /? , CKj <-> (3j , a <-> /? , a' <-> /?' 

so that after the elimination of the remaining parameters as a function of a, pj and p 
(j = 4, 5, • • • , N — 1), P^ m exhibits an ct-/3 symmetry 

iJH«,APi,p) = i3H/Wi,p) • (2-34) 

Similarly by the inversion of all variables Wi,Zi,U, r{ we have: 

IT (a, /?, Pi , p) = JJT(-2 - a - - 2p(s - 1) - p N - P, ftp,-, p) (2.35) 

where P = X] ^L/ Pi • Further information about IT can be obtained in the limit a — > oo 
(or /5 — > oo), by using a technique applied by Dotsenko and Fatteev 16 in the case of contour 
integrals, we found 

j-nm ^ a 2/3+2p(s-N-n+3)+2P-2m ^ 36) 

where we have used the kinematics: fci, k 2 , ■ ■ • , fcw-i > «o , k^ < &o and assumed cto < 
To eliminate most of the parameters as a function of a, f3,pj and p we use (2.4), (2.6) and 
momentum conservation. After such elimination the symmetry (2.25) becomes: 

IT (a, P,p j ,p) = I% n (m-l-P-a-P + p(N + n-l-s),P,p j ,p) (2.37) 

Using Stirling's formula, it is not difficult to check that the following Ansatz is consistent 
with (2.34), (2.36) and (2.37): 

AT = f^ m {p,p j )A{p-a)A(p-p)A(l-m+P+a+P+p( 8 +2-N-n)) 

AT = fTiP.Pj) fl A - h*)\ . (2 ' 38) 

j=i V / 
Now we can fix fTiPiPj) by using the 3-point function A§ m - 

8 N ~ 3 

A n N m (ki,k 2 ,k j ^Q,k N ) = (-Tr) N - 3 — Ar(ki,k 2 ,k N ) , 3<j<N-l .(2.39) 
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/aJV-3 \ N i 

(2.40) 



Now using the result for A^ m we get: 

)AT-3 \ ^ 

w ^- 3 )[A(-p)] 

i 

We finally return to (2.38) and obtain 

AT = (s + N - 3)(s + AT - 4) • • • (s + 1) [fuA(-p)] s 

N 1 

[- 7 rA(p- 1 )] m [-7rA(p)rn(- 7r ) A (2^-^)) ' ^ 

3 = 1 

therefore, redefining the screening operators, T^. and \i as before, we have: 

xN-3 

AT = ~Q^P S+N ~ 3 ( 2 - 42 ) 

which is a remarkable result; however, it is valid only in the kinematic region already 
mentioned. In order to extend for general k, we can use the same technique as used in 
[20]. Notice that the amplitude (2.41) factorizes as in the case without screening charges. 



3- Supersymmetric Correlators 

In a recent paper 19 we have calculated the 3- and 4-point NS correlations functions 
using DHK formulation 17 of super Liouville theory coupled to superconformal matter on 
the sphere without screening charges. The total action S is given by the sum of the super 
Liouville action Ssl and the matter piece Sm, 

Ssl = [ d 2 zE (\b a <S> S LD a ®SL - QY$ SL - 4ine a +* s A , 

I 3 \ (3 - 1} 

S M = J d 2 zE(^D^ M D a ^ M + 2ia Y§ M ) , 

where &m are super Liouville and matter superfields respectively. The central charge 
of the matter sector isc= §c, (c = 1 — 8o,q). Analogously to the bosonic case the 
parameters Q and a± are given by (compare with (2.2)) 

Q = 2^Jl + a 2 , a± = -| ± X - ^Q 2 - 4 = -| ± \a \ , «+«_ = 1 . (3.2) 

We shall call \I/ ns the gravitationally dressed primary superfields, whose form is given by 
V N s(zi,ki) = d 2 zEe ik ^ M ^e l3 ^ s ^ z \ where 

P(k) = -Q + \k-a \ . (3.3) 
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The calculation of the three-point function of the primary superfield S&ns, involves 
the expression: 

3 \ 3 

f[ J * N s{*iMV) = J [V^ SL ][V^ M ]l[^ NS (z t ,h)e- s . (3.4) 



Vi=l / i=l 



We closely follow the method already used in the bosonic case. After integrating over 
the bosonic zero modes we get 

Yl j *ats(z*, h^j = 2tvS k i - 2 «o j A 3 (h, fc 2 , fc 3 ) , 

(3°5) 

where (• • -)o denotes again the expectation value evaluated in the free theory (// = 0) and 
we have absorbed the factor [«+( — 7r/2) 3 ] -1 into the normalization of the path integral, 
the parameter s is defined as in the bosonic case (see (2.6)). 

For s non-negative integer, after fixing the SL2 gauge, 30 z 1 = , = 1 , Z3 = 00 , 62 = 
#3 = , 9\ = 9, in components ($sx = 4> + 9ifj + 9ip) (the integral above is the supersym- 
metric generalization of (B.9) of Ref.[16]) we have 

A(k 1 MM) = n-s)(^?( i ^ypi 

Z 7T 



s s 

X 

i=l i=l i<j 



« s s s 

j n n \zi\- 2a+pi \i-zi\- 2a+i32 n i^-^r 2a + ^(o)#(^) • • • w* a )> . 



(3.6) 

Observe that this is non- vanishing only for s odd (s = 21 + 1). One may evaluate (if) ■ ■ ■ ip)o 
and {ip • ■ -ip)o independently, since the rest of the integrand is symmetric, one may write 
the result in a simple form by relabelling coordinates: 

A,(k 1 MM) = n-s)(^f\C-^y^(-i)^s\\ 

s-l 

x j f[d 2 Zl n - n 1* - ^-i 4p n 1^-1 - ^r 2 k.r 2 

i=l i=l i<j i=l 

Redefining the variables as z s = w , -221- 1 = Ci an d z 2i = Vi we have: 

3 1 / 2 \ s 

^1,^2,^3) = -^ r(-s)r(s + l)— ) I z (a,/3 ;/ 9) , (3.7) 
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where 



I l (a,p;p) = ^-a 2 I d 2 w\[d 2 C i( i\ 

i=i 



\w 



2a-2 



i=i 



X 



nioi 2a ii-oi 2/3 ii 



Vi 



20 



i=l 



n & /' n & o " 

i,3 i<3 



Vi 



m\ 



i=l 



(3-8) 

and ct,/3,p are defined as before. In ref.[19] we calculated / in detail by using the sym- 
metries I l (a, (3; p) = I l (f3, a; p), I l (a, (3; p) = — a — (3 — 4/p, (3; p) and looking at its 
large a behavior we obtained: 



I l (a,f3;p) = - 



TV 



21+1 



2 2l 



21+1 I 



nA(2tp)lWi + (2t + l)p) 

i=l i=l ^ ' 



x Y[ A(l + a + 2ip)A(l + (3 + 2ip)A(-a - (3 + (2i - 41) p) 



i=0 
I 



x f[ A(i + a + (2i - l)p)A(i +(3+(2i-l)p)A(-±-a-(3+(2i-4l-l)p) 

i=l 

(3.9) 

We can choose, k\, ks > ao, &2 < «o- We proceed now as in the bosonic case, obtaining 
for the parameter /?, 

p(l - s) (a > 0) 
-| - ps (a < 0) . 



P 



(3.10) 



Now we are ready to write down the amplitude. For oto < we have the non-trivial 
amplitude: 



-Z7T, 



^(fcl,*2,*3) = (-^-) 3 



^ A Z^ 1 



- P 



1 s 



2 2 



-- A(l + a-(s- l)p) A - - a + p 



1 



n(-=)Afi[i +/ ?-fc 



(3.11) 



In the case c = 1, we obtain for the external legs, renormalization factors of the 
form A(l — \ki\), which should be compared to the bosonic case (2.30); it permits as well 
comparison to super matrix model as well, whenever those are available. 

By redefining the cosmological constant and the primary superfield ^ns 



Mi-/>) 



p , ^Ns(kj) 



1 



(-^)A(I[l + ^-fcJ]) 



we get 



As(ki, k 2 , k 3 ) = p 5 
13 



^Ns(kj) , (3.12) 



(3.13) 



As in the bosonic case we have a remarkably simple result. The only differences with 
respect to the bosonic case are in the details of the renormalization factors. Compare 
(3.12) with (2.11). Note that the singular point at the renormalization of the cosmological 
constant is p = — 1 in the bosonic case, which corresponds to c = 1, and p = — \ in the 
supersymmetric case, corresponding toc=lorc = 3/2. 

We shall now generalize the above result to the case which includes screening charges 
in the supermatter sector. We consider n charges e ld+ ® M and m charges e* d -* M , where 
d± are solutions of the equation ^d(d — 2ct ) = \- After integrating over the matter and 
Liouville zero modes we get 



i=l i=l i=l / 

= 2ivd ^2 k i + nd + + md ~ ~ 2a ^j ^3 m (^l 5 fc 2, fa) 

In „ ,2, m p ,2 

Ar(fa,fa,fa) = n-s)(^)H^) s (U J ^e^-^n J ^e"-*"<"> 



U=l i=l 



x 



J J^d 2 z i e iki ^ M ^e l3 ^ s ^ (J d 2 ze a +* SL ^ J , (3.14) 



Integrating over the Grasmann variables and fixing the SL(2) symmetry as before 
(zi = , Z2 = 1 , Z3 = oo , 6*i = 9 , 6*2 = 6*3 = 0) we obtain (using d + d- = —a+a- = —1) 



A nm {faM,fa) =T(-s 



tt\ 3 (ipa\\ s {-d\) n {-d 2 _y 



2 J \ 7r J n\ ml 

n p n 

xj[ dWur^i-ui-^niu-t, 

i=l J i<j 

m „ m n m 

x II / d 2 r t \r t \- 2d - kl \l ~ r t \- 2d ~ k2 U \ r * ~ 1111 ' r i 

xf[[ d 2 z i \z i \- 2 ^\l-z i \- 2a ^f[\z i -z j \- 2 < 

i=l i<j 
j n m s 

x / (^v(o) - kiim) n^(*o n^(^) n^(^) ) (3.15) 

\ i=i i=i i=i 



Since the vacuum expectation value of an odd number of ipip (or ) operators is zero 
we have only two non-trivial cases: in the first case n + m = odd, s = even and in the 
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second one n + m = even, s = odd. Thus we have (see also [21] for comparison) 



Ar(h,k 2 ,k 3 ) =T(-s 



where 



-ttV fi^lY (-d 2 ,) n (-d 2 _) 



x 



2 J \ 7T J n\ ml 
7™(a, (3; p) x /?; p) , n + m = even , s = odd 

J^P(a, (3;p) x J^(a, /?; p) , n + m = odd , s = even 



lb n lb 

i=i J i<j 



Tfh n TTL iv 1 1 o 

n / d'nin^'ii - r^' n \n - ^ nn ' 

i=l i<j i=l j=l 



m 



(3.16) 



\i=l i=l I 



/s s s I s \ 

n ^ n i^i 2 «ii - z % \^ n \* - ^ ( n^^o ) ( 3 - is ) 



with p = d+,p' = d 2 _ and a, (3, p, a, a' , (3, (3' defined as before. Note that I™ is the 
supersymmetric generalization of (B.10) of the second reference in [16]. The integral J™ 
differs from I^ 1 by the introduction of a factor ££(0) and Jq can be obtained from Iq by 
dropping ij)ij)(0). Henceforth we assume, for simplicity, n + m = even, s = odd. We will 
work out explicitly only the case n, m even. However, the final result for the amplitude 
does not depend on which case we choose. In ref. [24] we have calculated I^ 1 for n and 
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m even and we get: 



r^m - n^^,, (-I)" ""• [a (I - f)]" [a (i - |)] m 

2 ± 

x ~jf[ A(l + d + zp)A(l + /3 + zp) A(m - d - /3 + p(z - n + 1)) 

71 

x JjA(i + d + (z-^)p)A(^+^ + (z-^)p)A(-^-d + m-/3 + p(z-n + i)) 

7— 1 

777. ^ 

x JJ A(l + d' - ^ + ^p')A(l - J + ^ + i^)A(^ - d' - p' + p'{i - m + 1)) 

i=0 

m 

x H A(i-|+d' + (z-iy)A(^-^+^ + (z-^)A(-^+ ^-d'-^+p'(z-m+i)) 

i=i 

(3.19) 

In the case where s = 21 + 1 the gravitational contribution to ^.^(/ci, fo, ^3), i-e, Iq is 
just the same as in the case without screening charges, thus from the last section we have 
the supersymmetric generalization of (B.9) of ref.[16]: 



S-l 3-1 

nA(2zp)nA(i + (2z + l)p) 

i=l i=0 
J 

2 

x Yl A(l + d + 2zp)A(l + (3 + 2zp)A(-d -(3 + 2p(i - s + 1)) 

z=0 
3-1 

x f[ A(i + d + (2z - l)p)A(^ + + (2i - l)p)A(-i - d - £ + p(2z - 2s + 1)) 

i=0 

(3.20) 

To obtain A^ m (ki, k 2 , k 3 ) (see (3.16)) we have to calculate I nm x I G . Using the same 
kinematics as in the case without screening charges, and considering that a < it is easy 



1 G 



_l TV' 

2^=1' 



A.--P 



16 



to deduce: 



a = a — 2p , a! = — 1 + 




P = ~-j-(n + s)p (3.21) 
P=-p-l = (n+l)p+^-l 

^' = ^U^V- 1) • 



Substituting in (3.19) and (3.20) and using (3.16) we obtain a very involved expression 
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(see also [26,21]): 

= r(- 8 ) 



n 



n 



*M 2(s-l) 



(2p) n - m (-)^ 



ra + in + s-1 7T 



s+n+m 



-p- 2mn sl 



X 



I ~ -n- . , n ip 



x J] A(-2ip)A(- + p(l - 2i)) J] A(-- 



i=l 



2 2 V4 2 



x 



X 



X 



II A G + Y + P^ + s -2^))riA(y + p (n + S + l-2z)) 

i=0 ^ ' i=l 

m t rn 

1 s p" 1 



i=0 



1 n, 

2 ~ 2 



2 2 



2 ('-T 



A(--p) 



3-1 S-1 

fl A(2zp) n A (i + (2i + l)p) 

i=l i=0 ^ ' 



s-1 
2 



X 



nA(i-^ + (2z-n- S )p)nA(- 



i=0 

n 
2 



(2i — 1 — n — s)p) 



i=i 



1 1 

x H A(- + a - p{2i + 1))A(- + ^- a -p(s-n-2 + 2i)) 
i=i 

i-i i-i 

2 1 2 1 

xJ]A(^ + -- a + (2z- S + n + 2)p) J] A(- + a + (2i - l)p) 



i=0 

71 

2 



'Y[ A(y - a - p(2i -n + s- 1))A(1 + a - 2zp) 



X 



i=l 

3-1 

2 



3-1 

2 



x H A(l + a + 2ip) Y[ A(^ - a + (2i - s + n + l)p) 



i=0 

rn 
2 



i=l 



X 



i=l 

rn 
2 



n - | - £(i + « - £))A(-£ - V(' " 1 " «)) 



-1 

2 2~ 



x 



IIA(- 



^-l-^(z-a-l))A(l-^-^(z- (m + 1) 
2 2 2 v 2^ v 2 2 v 2 



+ «)) 



(3.22) 



In order to obtain a simple expression for the amplitude we have to combine in each 
term the matter and the gravitational parts as in the bosonic case. The calculation is more 
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complicate now, but we finally get 

3 



ATikiMM = (-| 



2 v 2 ' > 



ITT . . 1 £> 1 

— A(- + - — 

2 v 2 4 



— A(- +p) 

2 v 2 p; 



xA p-a + - A - 



1 n + s mp 1 



A(l - — + a + (s-n-l)p) 



/i 



1 



x 



n 



2 

Z7T 



Z7T 



^A ^ + 



1 , P 



2 4 

2\ 



%7T 



(3.23) 



Therefore after redefining the cosmo logical constant, the NS operators and the screening 
charges 



id + <S> M {ti) 



id_<t> M {ti) 



in 



1 



-i -l 



in 



^A ^ + 



1 , P 



id + $ M (ti) 



e 

-l 



NS 



P 



in 



1 1 



_A - + -(#- fc; 



n -1 



JV5 



-1 



/i 



we obtain the very simple result: 



(3.24a) 
(3.246) 
(3.24c) 
(3.24a 1 ) 

(3.25) 



In view of the complexity of (3.22), the simplicity of the result is remarkable. 

As in the bosonic case we can calculate ratios of correlation functions either using 
(3.23) or (3.25). We obtain in a generic kinematic region: 



R = 



ULl(^Ns(k rj r')^Ns(k rj r')) 



R = (2a+) 



s EUi \nd- +rjd+\ 
(2p-l)(2p+l) 



(3.26) 
(3.27) 



Compare with (2.27). 

The above result agrees with other results 26 ' 21 simultaneously obtained in the litera- 
ture. Although the the continuations to non integer values of s used in [26] and [21] are 
not the same and, in principle, do not correspond to the procedure used here, the physical 
results seem to be independent of such details. 
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We now show that it is possible to obtain a simple result for the most general case of a 
iV-point amplitude with an arbitrary number of screening charges (A 7 }] 71 ). In that general 
case to compute the amplitudes we have to calculate the following integral 

AT(ki, ■■-,k N ) = ' (j[ J £?2 4e 'fc*"( a O+A*«(iO 

x f[J d 2 Ue id ^ M ^ f[ J d 2 r je id -* M ^ f[ J d\e a +* SL ^ 



where s = —-^(YliLi A + Q) an d Si=i hi + n( ^+ + fnd- = 2ct - After fixing the SL 2 
symmetry as before and integrating over the Grassmann variables the amplitude becomes 



^ r= r(- S )(-|) 3 (^)V4r(-a 



\ 7T J 

ii / ^n / ^n / ^nrfVkii- 2a+ ^ii-^r 2Q+/32 

7=4^ £=]V £=lV i=l 



N 

X 

3- 

s N 

x n k - w j\~ 2a% n n k - ^r 2a+ft 

i<j i=lj=4 
N N 

j=4 j<;=4 

n n n m 

x n i^i 2fc2d+ ii - ^i 2fc2d + n i** - *ii 2<p+ nn >■ ^- 2 

j=l i<j i=lj = l 

m m n N m N 

x n n 2m -ii - r,i 2fc2d - n in - ^i 2d - nn^- ~^\ 2d+kj nn ' ^i 2d - fcj 

j=l i<j i=l j=4 i=lj=4 

/ JV n m s \ 

x /( / 3 1 2 ^(o)-fc 2 ^(o)) n(^(%)-^ 2 ee(%)) n^^)n^^)n^(^)/ 3 - 29 ) 

\ j=4 j=l 1=1 1=1 I Q 

Now we have several terms which give non-trivial amplitudes, in the following we assume 
m + n and N + s even; thus we have: 



3/' 2 \ s I N \ 

N „ n „ m „ s 

x« 2 n / / ^n / ^jj^vKi^ii-wii^ 

j=4^ i=lV i=]V i=l 

s N N N 

xnK-^i^nnK-^i^ni^i 20 '! 1 -^! 2 ^ n i^-^ 

i<j i=lj=4 j=4 J<i=4 
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X 



n i^i 23 ! i - **i 2 ^ n i** - *^i 2 ^ n ii i** - 

i=l i<j i=lj = l 

m m n N 

x n i^i 23 '! 1 - ^ n ^ - ^ * n n - ^ n n \* - * 

i=lj=4 



m N 
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i=l j=4 



Is N 

x (n«(*on«( r *)) #(o)n#wn^fe) 

\ J = 1 J= 4 



(3.30) 



Vi=l 



i=l 



The definitions of the kinematics parameters are defined as in the bosonic case 

We shall use the kinematics: k±, &2, • • • , fejv-i > «o , /cat < cto < in order to eliminate 

all parameters in terms of a, /?, p and pj (4 < j < iV — 1). 
The symmetries: 



A N ™(a,p,p,p 1 ,p 2 , • • • ,Pjv-i) = -4aT(/3, a,p,Pi,P2, • • • ,PAT-l) 



(3.31) 



^ m (a,A p,Pi • • -pat-i) =A n N m {-^l3+ y -—^--P+ p(N+n- 8 -l), 0, p, Pl >- -piv-i) 



AT-l 

with P = X) Pj and the large-a behaviour: A% m (a -> oo) ~ a i-™+2/3+2 P ( s -7V-n+3)+2P 

motivate us to write down the ansatz: 

•AST = /HP, Pi, • • • , Pn-i) A(i+p-a) A(i+p-$ A(l - y +P+a+/3+p(2+ S -n-A)) (3.32) 

Taking the limit fcj — > (3 < z < A" — 1), which implies pj — > 2p (4 < j < A" — 1), we can 
determine /# m (p,pi, ■ ■ -pjv-i) using: 

2 J a, 



(3.33) 



and the result for A^ m (see (3.23)). We arrive at the result 



Anm 



(s + N - 3)(s + N - 4) ■ ■ ■ (s + 1) 



/"A(--p) 



x 



"1" A / 1 \ 

A(- + p) 

2 v 2 p; 



--A(i + ^) 
2 V 2 4 ; 



i N 

n 

i=i 



AT 



~) A(i(l + A 2 -^ 2 ))(3.34) 



Redefining ^>ns > A 4 an d the screening charges we have our final result: 



/inm 



AT-3 



fj, 



s+iV-3 



(3.35) 



which has the same functional form as the bosonic amplitude (2.42). As in that case the 
above result is correct only in the kinematic region used to calculate it and has to be 
continued outside this region. 
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4- Conclusion 



In the first part of this paper we have generalized previous results 9 for the Appoint 
tachyon correlator in Liouville theory coupled to conformal matter (on the sphere) with 
c < 1 to the case which includes screening charges in the matter sector. The results 
might be useful in understanding the issue of fusion rules in the calculation of the 3-point 
correlator (see discussion in [20]). 

In the second part we have obtained the A-point NS-correlators in super Liouville 
coupled to c < 1 matter (also on the sphere), including screening charges thus generalizing 
the results of [21,26] to the limit case c = 1, (N > 3), and the results of [20], obtained for 
s = 0, to any value of s. 

We have obtained, explicit formulae for the corresponding 2D-integrals involved and 
the final form of the amplitudes factorizes in the iV-external legs factors, confirming the 
results obtained in [20] (for s = 0) through a detailed analysis of the pole structure of the 
integrals. 

In our calculations it was possible to see the singularity in the renormalization of the 
cosmological constant (p — > /2/A(| — p)) at the point c = 1 (p = —1/2). This is similar 
to the bosonic case where p — > p/A(—p) and the c = 1 , (p = —1) point is also singular. 
The final (renormalized) A-point amplitude has the same form of the bosonic one. The 
similarity with the bosonic case has been found before in the discrete approach [20,25]. 

Finally we should stress that our results must be continued to other kinematic regions, 
this is likely to be very similar to the bosonic case (without s.c.) worked out in [20]. 

As a further development it would be interesting to carry out analogous calculations 
in the case with N = 2 supersymmetry and see, among other aspects, whether the barrier 
at c = 3 indeed disappears. The inclusion of the Ramond sector (for s ^ 0) is also of 
interest. 
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